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We show how two level atoms can be used to determine the local time dependent spectrum. 
The method is applied to a one dimensional cavity. The spectrum obtained is compared with the 
mode spectrum determined using spatially filtered second order correlation functions. The spectra 
obtained using two level atoms give identical results with the mode spectrum. One benefit of the 
method is that only one time averages are needed. It is also more closely related to a realistic 
measurement scheme than any other definition of a time dependent spectrum. 



I. INTRODUCTION 

It has turned out to be difRcult to give a universally ac- 
cepted definition of the time dependent spectrum. There 
are several definitions jl]-H] which have not been popu- 
lar in optics. A definition, which is connected to a re- 
alistic spectrum measurement was given by Eberly and 
Wodkiewicz Their, 'physical spectrum', has become 
a kind of canonical definition for a time dependent spec- 
trum in quantum optics. All definitions mentioned above 
are based on integrals, typically Fourier transforms, over 
different two time averages of classical stochastic vari- 
ables or quantum mechanical operators. 

In our earlier paper a totally new approach was 
introduced. Instead of calculating two time averages the 
radiation is guided to a group of two level atoms. All 
the atoms have equal, very small decay constants, but 
the resonance frequencies are all different. The spectrum 
measured by the atoms can be read from the excitation 
probabilities of the atoms. We have studied the spectrum 
of resonance fluorescence radiation from a laser driven 
three level atom ||]. The method was shown to give 
exactly the same result as the 'physical spectrum' de- 
fined by Eberly and Wodkiewicz. In order to calculate 
the two time averages needed to determine the 'physical 
spectrum', the Quantum Regression Theorem (QRT) ||] 
must be used. It has been shown that there are systems 
where QRT cannot be used i.e. the two time averages it 
gives are incorrect . This is the case if the interaction 
between the matter, which emits the radiation, and the 
radiation field is strong. Interaction can be strong e.g. 
in microcavities. In the method of analyzer atoms, two 
time averages are not needed. The method can be used 
also in situations where the QRT is known to give incor- 
rect results. It can be said that it is closer to a realistic 
spectrum measurement than any other definition. 

In our earlier paper, the three level atom was coupled 
to the analyzer atoms using cascaded master equations 
Thus the quantum mechanical state of the field 
was not made available. In this paper we connect the 
method of analyzer atoms to our cavity QED simulations 
in one dimension. In these simulations there are two level 
atoms inside a one dimensional cavity. The field is quan- 



tized using canonical quantization. The state vector is re- 
stricted to have only one excitation i.e. the field strength 
is very low. The model was introduced by V. Buzek et. 
al. |10| , [11[ and has been used in several different studies 
in one | [12[ and two dimensions | ]l3t . The quantum me- 
chanical state of the field is available so the spectrum of 
the radiation can be read directly from the excitations 
of the modes. The spatial dependence of the spectrum 
is obtained by calculating the second order correlation 
function g{ri,r2) and using appropriate filter functions. 
The typical situation in the simulations is that the ini- 
tial field intensity is split into left and right propagating 
parts using two level atoms. The spectra of both parts is 
determined using both analyzer atoms and filtered cor- 
relation functions. In all cases the two different methods 
give identical results. 

In Sec. II we specify the model used in the simula- 
tions. Spatial correlation functions and their relation to 
the spectrum is explained in Sec. III. In Sec. IV we 
present the analyzer atom method for spectrum mea- 
surements. In Sec. IV we show several results of our 
simulations and finally in Sec. VI we present our conclu- 



II. HAMILTONIANS AND STATE VECTORS 

A. Hamiltonians 

In our simulations the mode functions in the one di- 
mensional cavity (0 < r < L) are 
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G'„(r) sin(fc„r), fc„ = n— , ?7 = 1,2,3... (1) 

All mode functions are zero at the edges of the cavity, 
r = and r = L. The electric and magnetic field opera- 
tors can be expanded using the mode functions (|^) 
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where a and are annihilation and creation operators. 
The operators satisfy the usual canonical commutation 
relations 
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For the energy density operator we get using the ex- 
pansions (0) 
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Integration of H{r) over the whole cavity and use of the 
orthogonal integral 
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sin(fc„r) sin(fc„/r)dr = — (5„™, rt, m > 1 (6) 



gives the familiar field Hamiltonian 
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At fixed positions inside the cavity there are Na two 
level atoms with resonance frequencies uoj and dipole con- 
stants Dj. The atomic Hamiltonian is the sum over all 
one atom Hamiltonians 
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where dz is a Pauli spin matrix. The radiation field 
and the atoms are coupled through dipole coupling. The 
dipole operator of the j:th atom is 
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For the interaction Hamiltonian we get 
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where E{rj) is the electric field operator (g) at the atomic 
position Tj . Using the expansion (|^) for the field operator 
and (0) for the dipole operator we get 
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The terms CT+aj and a^aj, which do not affect the time 
evolution significantly, have been neglected. This ap- 
proximation is called the Rotating Wave approximation 
(RWA). The Hamiltonian ( pr] ) has the familiar Jaynes- 
Cummings form. The total Hamiltonian is the sum of 
field, atomic and interaction Hamiltonians: 



H — Hp + Ha + Hj . 



(12) 




B. State vectors 



The general state vector is restricted to have only one 
excitation. The most general state vector has the form 

Na 

k' \'^Y[ \0)j 
J = l 
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In the first sum, one of the field basis vectors have one 
excitation and in the second one, the excitation is in the 
atomic basis functions. If the initial state vector has 
only one excitation it does not obtain more excitations 
because the RWA was used in the interaction Hamilto- 
nian. The state vector at all times is given by equation 

(HI)- 

The state vector used as the initial state in our simu- 
lations is of the form 
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The mode distribution |cfe|^ is a Gaussian centered at 
k = ko with a variance cr^.. The phase factor e"*'^''" is 
important. It guarantees that the energy density distri- 
bution is also a Gaussian centered at tq. 



III. CORRELATION FUNCTIONS AND THE 
MODE SPECTRUM 



A. Correlation functions 



The exact quantum mechanical state of the field is de- 
termined uniquely if all correlation functions of the oper- 
ators E{r) and B(r) are known. In our simulations, be- 
cause of the special form of the state vector (|l^ ) the sec- 
ond order correlation function determines the field state 
uniquely. In the following we denote ~ fiQ — 1. Us- 
ing the expansions (|^) we get for the normally ordered 
correlation function 



2 



: E{ri)E{r2) : 

= ^ ^ j^uJnUJm sin(/c„ri) sin(fc„ir2)(5'(i)|at a„ + a]„a„|^'(i)) 



(15) 



n— 1 m— 1 

oo oo 



— V^n^m sin(A;nri) s\n{kmr2){c^{t)cm{t) + cj;,(i)cn(t)) 



n— 1 m— 1 

T*(ri,i)T(r2,t)+T(ri,i)T*(r2,t) 



where 



'T'C^ ^) = X! V ~r sin(A:„r)c„(t) 



(16) 



A similar calculation gives 



{^{t)\ : B(ri)S(r2) : \^{t)) (17) 
= r*(ri,t)r(r2,i)-r*(r2,<)r(ri,t). 

It is possible to invert the formula (|6|) and get the 
mode coefficients in Eq. ([l^). Multiplying both sides 
by sin(fcmr) and using the orthogonality integral (^) we 
get 
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So if T(r, t) is known the mode coefficients can be cal- 
culated. It is also possible to calculate the coefficients 
c„ from the correlation functions ( p^ ) and (pT[). Let us 
define 

W^(ri,r2) 

= {^{t)\ (: ^(ri)i?(r2) : + : B{r,)B{r2) :) |vl'(t)> 
= T*{n)T{r2) (19) 
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Integrating both sides by J sin(fcir) sin(fcpr)dr we get 



to get a spectrum of radiation in some part of the cavity 
only we use filtered correlation functions. The correlation 
function ( [l9| ) is replaced by a filtered correlation function 



WF{ri,r2;ro) = {ri,r2,rQ)W {ri,r2), 



(21) 



where the filter function is a real- valued window function 
centered at rg. The method is similar to the Windowed 
Fourier Transform (WFT), Ref. which is used to 

determine the time dependent frequency distribution of a 
time dependent signal. The mode spectrum is calculated 
using a filtered correlation function (|2l|) in equation ([2^). 
The spectrum obtained is the mode spectrum of the ra- 
diation in the part of the cavity where the filter function 
is nonzero. 

A typical example of the filter function of one variable 
is a Gaussian filter 
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where tr^ is the variance of the window. Two dimensional 
filter is a product of two one dimensional ones. Another 
filter used in our simulations is a constant filter which is 
unity in some region and zero everywhere else 
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The general behavior of the spectrum with all filter func- 
tions is such that the broader the filter function, the 
smaller details can be seen in the spectrum. 
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This gives equations for each Cp . For p = 1 we get | ci ^ . 
We can choose ci to be real and determine ci. The rest 
of the coefficients are now determined uniquely. Normal- 
ization can be used as a check of the calculation. 



B. A local mode spectrum 

The mode spectrum |cfep gives the time dependent 
spectrum of the radiation in the whole cavity. In order 



IV. TRADITIONAL DEFINITIONS OF THE 
SPECTRUM AND ANALYZER ATOMS 

A. Traditional definitions of the spectrum 

In the previous section the spectrum was calculated us- 
ing ^spatial correlation functions at a specific time value\ 
The traditional method to get the spectrum is to calcu- 
late the Fourier transform of two time averages of field 
operators, i.e. use Hime correlation functions at a spe- 
cific point. For a stationary field the spectrum can be 
defined as 



S{uj) = / dT{E{t)E{t + T))e" 



(24) 
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This definition is valid only for stationary fields. There Eberly and Wodkiewicz They define a 'physical spec- 
are many generalizations of a Fourier spectrum for non- trum' to be a double Fourier transform of two time aver- 
stationary fields [0-^. A method which takes also ages multiplied by filter functions 
the measurement scheme into account was presented by 

oo oo 




The frequency filters are placed in front of the pho- 
todetector. They allow only the radiation with a certain 
frequency to pass to the detector. The spectrum is ob- 
tained from the relative intensities measured by the de- 
tector when the filters are tuned to different frequencies. 
The spectrum obtained is time dependent and the defini- 
tion incorporates the measurement scheme. Our method 
to include time dependence to the spectrum is analogous 
to the method shown in the previous section. The defini- 
tion ( p5| ) has 'time filtered^ correlation functions whereas 
in the last section \spatially filtered^ correlation functions 
were used. 

In many problems of quantum optics the radiation field 
is traced out from the Hilbert space. As a result the time 
evolution of the system is described by a master equation, 
which does not have information about the field degrees 
of freedom. Using this approach it is not possible to 
calculate the spatial correlation functions; then the defi- 
nitions with time correlation functions like (|2^) and ( |25| ) 
are the only useful definitions. 



B. Analyzer atoms 

In our earlier paper |^ , a totally different approach to 
the time dependent spectrum was introduced. The radi- 
ation from the system of interest is guided to a group of 
N two level atoms which all have the same very small de- 
cay constant F. The resonance frequencies of the atoms 
are all different 

A A ^max ^min -i r» j\r 

UJ, lAU) = , n=l,2...iv. 

N-1 

(26) 

Initially all atoms are in the ground state. The incoming 
radiation excites the atoms. Because the decay constants 
are equal and very small, the excitation probabilities of 
the atoms are directly proportional to the intensity of the 
incoming radiation at the resonance frequency. The nor- 
malized excitation as a function of resonance frequency 
of the atoms gives the time dependent spectrum of the 
radiation. A detailed description of the method and a 
comparison with a spectrum (^5|) calculated using two 
time averages can be found in Ref. In these simu- 
lations cascaded master equations were used to describe 
a three level atom and a group of analyzer atoms. The 



two spectra were shown to give identical results for the 
resonance fluorescence radiation of a three level atom. 

The analyzer atom method can be used in our cavity 
simulations. Now the atoms with frequencies given by 
Eq. (^6|) and very small decay constants are located at 
some point inside the cavity. As in the case of master 
equations, the excitations should give the local time de- 
pendent spectrum of the radiation which passes them. In 
the next section the spectrum measured using detector 
atoms is compared with the mode spectrum calculated 
using filtered spatial correlation functions. 

V. SIMULATIONS 

In this section we show results of spectrum measure- 
ments in three different cases. In all simulations there 
is one atom or several atoms at the center of the cav- 
ity. The initial field state is localized on the left and is 
moving to the right. The atom at the center splits the 
photon into left and right propagating parts. The spectra 
of these two parts are measured separately using analyzer 
atoms and filtered correlation functions. In the first sim- 
ulation the initial field state is the Gaussian (^4|) and 
there is one atom at the center. In the second simulation 
there are three atoms at the center. In the last simula- 
tion there is only one center atom and the field state is 
a superposition of several random Gaussian states. The 
units in our simulations are chosen is such a way that 
c = eo = fJ.0 = h = 1. 

A. A Gaussian photon and one atom 

In the first simulation a Gaussian photon propa- 
gates towards the center of the cavity. The length of the 
cavity is i = 27r. The mode distribution is a Gaussian 
(Q), centered at fco = 100.0 with the width Tph = 47r. 
The intensity profile is centered at rg = 2.0. The initial 
intensity is shown in Fig. |l|. 
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FIG. 1. The energy density of the photon at t = 0.0 and 
t — 3.8. The initial energy density is a Gaussian (^4|) with 
parameters ro = 2.0, ko — 100.0 and Ft = An. At t = 3.8 
there is one peak propagating to the left and two to the right. 
The initial intensity is split as a results of the interaction 
with the center atom, marked by a cross. The circles show 
the positions of the analyzer atoms. The length of the cavity 
is L = 27r and the number of modes is Nmode = 400. 

In the center at r = , there is a two level atom vifhich 
is exactly on resonance with the photon cjo = 100.0. The 
decay constant of the atom is F = tt, so the linewidth of 
the atom is narrower than the width of the mode distri- 
bution of the field. When the photon reaches the atom, 
the atom gets a nonzero excitation. Later the atom emits 
energy back to the field modes via decay to the ground 
state. The population of the excited state as a function 
of time is shown in Fig. |^. 




As a result of the interaction, part of the energy is re- 
flected to the left and part is able to pass the atom. The 
intensity &tt — 3.8, after the interaction, is shown in Fig. 
|l|. The intensity profile on the right has two peaks. The 
first peak has propagated to the right without interaction 
and the second is the result of the atomic decay. 

In addition to the center atom there are 600 atoms 
which detect the spectrum of the radiation which passes 
them as was explained in section IV. At r = 1.8 i.e. left 
from the center atom there are 200 atoms which mea- 
sure the spectrum of the reflected radiation. The dipole 
couplings of the atoms are time dependent in such a way 
that they are zero when the initial Gaussian photon prop- 
agates to the right. Kit = 1.5 the atomic dipoles get the 
values determined by and start to detect the spec- 
trum of the radiation. On the right at r = -I- 1.0 there 
are two groups of analyzer atoms. The first group (200 
atoms) measures the spectrum all the time. After the 
first peak on the right has passed them, its spectrum 
can be read from the excitations. After the second peak 
has passed them the atoms give the spectrum of all ra- 
diation right from the atom i.e. spectrum of the first 
and second peaks. Another group at the same position 
r = L/2.0 -f- 1.0 detects the spectrum of the second peak 
on the right. The dipole constants of these atoms get 
nonzero values ( p6| ) immediately after the first peak has 
passed them. 

The spectra of the total intensity on the right and left 
measured using the atoms are shown in Fig. 0. 




FIG. 2. The excitation of the center atom. The photon 
excites the atom which decays exponentially to the ground 
state. The decay constant of the atom is F = tt. 



FIG. 3. The normalized spectra of the radiation on the 
left and right measured by the analyzer atoms. Fig. ^. The 
resonant radiation has been reflected to the left and the 
off-resonant radiation has propagated to the right. The spec- 
trum has been read from the atoms after all radiation has 
propagated through them. Both spectra have been normal- 
ized in such a way that the area under the curves is unity. 
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On the right the spectrum has a two peak structure. 
It is interesting that on the right there is no intensity at 
the resonance frequency of the atom. The center atom 
has been able to reflect the resonant radiation to the left 
and only off-resonant radiation passes the atom. The fact 
that resonant radiation is missing is interesting because 
the second peak of the intensity profile is the result of 
the atomic decay. 

As was explained earlier, we also measured the spectra 
of the two peaks on the right separately. The result is 
shown in Fig. ^. 



First peak 
Second peak 




FIG. 4. The measured spectra of the two peaks on the right 
(Fig. |l|) separately. The spectrum of the first peak is broad 
because the intensity profile is narrow. The second peak has 
a typical Lorentzian spectrum of free decay. When the spec- 
trum of both peaks is measured, as in Fig. |^, the second 
peak creates an stimulated decay type of effect to the detec- 
tor atoms and the measured spectrum has a dip in the middle. 
The normalization of the spectra is the same as in Fig. ^ 

The spectrum of the first peak is quite close to the 
initial spectrum. The spectrum of the second peak is 
narrower and is similar to the spectrum of the sponta- 
neous decay. The fact that the spectrum of the second 
peak is narrower is understandable because the intensity 
profile is wider and the interaction time with the detector 
atoms longer. The combined spectrum of the two peaks 
has a dip in the middle. Fig. ^, whereas the spectra of 
the peaks measured separately both have a single peak 
structure. This shows that the spectrum is not additive. 
The second peak causes a stimulated decay type of effect 
and the excitations of the detector atoms, which are on 
resonance with the radiation, get smaller. 

Next we calculate the mode spectrum using the filtered 
correlation functions as was explained in Sec. III. The 
normally ordered correlation function (: E{ri)E{r2) :) at 
t = 3.8 is shown in Fig. || (a). The diagonal elements 
Ti = T2 give the energy density of the radiation. Off- 



diagonal elements show the coherence of the radiation. 
In order to calculate the spectrum of the radiation in 
the left part of the cavity, a filtered correlation function 
is used. The part ri,r2 > ^ is replaced by zero. The 
filtered correlation function is shown in Fig. ^ (b). 
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FIG. 5. The correlation function {*! : E{ri)E{r2) : |*> of 
the radiation a.t t = 3.8 (a). The diagonal elements give the 
energy density of the radiation, Fig. |l|. The lower figure (b) 
shows the filtered correlation function used to determine the 
mode spectrum on the right. 

Similarly in order to get the spectrum of the radia- 
tion on the right, the replacement (: E{ri)E{r2) :) = 
when ri,r2 < is used. The mode spectrum is calcu- 
lated using the formula (^0|) with a filtered correlation 
function. The mode spectra obtained are identical to the 
spectra measured using the analyzer atoms. Fig. ||. We 
calculated also the mode spectra of the two peaks on the 
right separately using the appropriate filtered correlation 
functions. Also these spectra were the same as measured 
by analyzer atoms. Fig. 0. 
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B. A Gaussian photon and three atoms 



C. A random photon and one atom 



Next we add two more atoms to the center with reso- 
nance frequencies and decay constants uo = 90.0, T — tt 
and ujQ ~ 110.0, F = 7r/4 respectively. The third atom 
is the same as earUer, luq — 100.0 and F = tt. The atom 
with resonance frequency = 110.0 has a smaller de- 
cay constant than the other two atoms. The width of 
the mode spectrum of the Gaussian initial photon is in- 
creased to Fp/i — 8tt. Because the spectrum is very broad, 
the energy density profile of the photon is narrow. We 
have carried out exactly the same spectrum measurement 
as in the previous simulation for the radiation which 
passes the atoms and for the reflected radiation. The 
results are shown in Fig. ^ 



In the two previous simulations the initial state of the 
photon has been of the Gaussian form ([T^). Next we 
choose the initial state to be rather exotic. The idea of 
this simulation is to test the method using a totally dif- 
ferent type of field excitation than in the previous sim- 
ulations. The initial photon is a superposition of ten 
random Gaussian states. The parameters kg, of, and 
in equation (|l^) are chosen randomly, in such a way that 
the initial energy density is on the left and propagating 
to the right. The frequency distribution is centered at 
Lo — 100.0. Again at the center there is one atom which 
splits the intensity into two parts. The initial intensity 
profile of the photon is shown in Fig. 0(a). 




FIG. 6. The measured spectra of the radiation on the left 
and right in the simulation when there are three atoms at 
the center of the cavity. The radiation on the left has three 
peaks at resonance frequencies. On the right the broad spec- 
trum has three holes. The width of the peak (or hole) at 
uJo = 110.0 is narrower than the two others. The length of 
the cavity is L = Stt and the number of modes Nmode ~ 1600. 
The normalization of the spectra is the same as in Fig. ^ 

On the left the spectrum has three peaks at the res- 
onance frequency. On the right the original Gaussian 
spectrum has three holes at the atomic frequencies. The 
width of the peak and the hole at = HO is narrower 
as was expected based on the atomic parameters. The 
result is qualitatively the same as in the simulation with 
one atom in the center. The resonant radiation is re- 
flected and the off-resonant radiation is able to pass the 
atoms. As in the previous simulation we calculated the 
two spectra using also the flltered correlation functions. 
The correlation function has a more complicated form 
than the one in Fig. |^. The result is again identical to 
the spectra given by the analyzer atoms. 
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FIG. 7. The intensity profile (a) and the mode spectrum 
(b) of the field state which is a sum of ten Gaussian distri- 
butions for the photon. The field is propagating to the right. 
Note that only a part of the cavity is shown. At the cen- 
ter, = ■§■, there is one atom with the resonance frequency 
LJo = 100.0. The circles on the left and right are analyzer 
atom which detect the spectrum. The initial mode spectrum, 
(b), has a six peak structure as some of the Gaussians are 
overlapping. The normalization of the spectrum is the same 
as in Fig. g 

Note that the length of the cavity is now L — 8tt. The 
fleld has a four peak structure and it propagates to the 
right towards the center atom (thick circle). The atom 
has the resonance frequency ujq = 100.0. Circles right 
and left from the center atom show the positions of the 
analyzer atoms which detect the spectrum. Figure ^(b) 
shows the mode spectrum of the initial photon. 

As in the earlier simulations the center atom splits the 
radiation into left and right propagating parts. The spec- 
tra of the two parts are shown in Fig. @. 
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frequency 

FIG. 8. The spectra on the left and right measured using 
analyzer atoms. The initial state of the field is a sum of ten 
random Gaussian states, Fig. ^. The spectrum has been read 
from the atoms after all radiation has propagated through 
them. The normalization of the spectrum is the same as in 
Fig. I 

Because there is not much intensity which would be 
at resonance with the center atom, the spectrum on the 
right is similar to the original mode spectrum. Only a 
small fraction of the intensity is reflected to the left. Both 
spectra have several peaks. We calculated the mode spec- 
trum on the left and right separately using filter functions 
as earlier. Again the mode spectrum gives exactly the 
same results as detected using analyzer atoms. 



VI. CONCLUSION 

We have used two level atoms to detect the time de- 
pendent and local spectrum of radiation in a one dimen- 
sional cavity. The spectrum is read from the excitation 
probabilities of atoms with very small decay constants. 
An alternative method to determine the spectrum in this 
paper has been to calculate the mode coefficients using a 
filtered second order correlation function. Different filter 
functions can be used to determine the spectrum at a 
specific part of the cavity. It might be possible to gener- 
alize the approach of filtered correlation functions using 
wavelet expansions ||-|6). InaU cases studied, the spec- 
tra determined using these two methods give identical 
results. In our earlier paper we showed that the analyzer 
atom spectrum gives the same result for a resonance fluo- 



rescence spectrum of a laser driven three level atom as the 
'physical spectrum' defined by Eberly and Wodkiewicz. 
This paper gives further proof that the method really 
works. 

One benefit of the method is that only one time aver- 
ages of quantum mechanical operators are needed. Thus 
it is closer to a realistic spectrum measurement than the 
usual definitions of time dependent spectrum, which typ- 
ically require two time averages. It can also be used in 
situations were the usual method to calculate two time 
averages in the Schrodinger picture, the quantum regres- 
sion theorem, is known to give incorrect results. 
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